Auction-based Resource Allocation in 
MillimeterWave Wireless Access Networks 

George Athanasiou, Member, IEEE, Pradeep Chathuranga 
Weeraddana, Member, IEEE, and Carlo Fischione, Member, IEEE, 
&) 

o 

O ' Abstract 

< 

f. The resource allocation problem of optimal assignment of the clients to the available access points 

in 60 GHz millimeterWave wireless access networks is investigated. The problem is posed as a multi- 
assignment optimization problem. The proposed solution method converts the initial problem to a 
minimum cost flow problem and allows to design an efficient algorithm by a combination of auction 

O ■ 

algorithms. The solution algorithm exploits the network optimization structure of the problem, and thus is 
much more powerful than computationally intensive general-purpose solvers. Theoretical and numerical 
results evince numerous properties, such as optimality, convergence, and scalability in comparison to 
existing approaches. 
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I. Introduction 
MillimetterWave (mmW) communications utilize the part of the electromagnetic spectrum 
between 30 and 300 GHz, which corresponds to wavelengths from 10 mm to 1 mm [QQ|. MmW 
^ wireless networks in the 60 GHz unlicensed band are considered one of the key technologies 

S-j ■ 

for enabling multi-gigabit wireless access (transmission rates up to 7 Gbps) and provisioning 
of QoS-sensitive applications. Multiple industry-led efforts and international organizations have 
emerged for the standardization . More than 5 GHz of continuous bandwidth is available in many 
countries worldwide, which makes 60 GHz systems particularly attractive for gigabit wireless 
applications such as gigabyte file transfer, wireless docking station, wireless gigabit ethernet, 
wireless gaming, and uncompressed high definition video transmission. Moreover, scenarios such 
as dense small-cells and mobile data offloading [2], which are nowadays strongly motivated by 
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the increased end-user connectivity requirements and mobile traffic, can be accommodated with 
the use of 60 GHz radio access technology. 

Resource allocation for wireless local area networks has been the focus of intense research. 
Several studies have analyzed the performance of the basic client association policy that IEEE 
802.11 standard defines, based on the received signal strength indicator (RSSI). These studies 
have showed that this basic association policy can lead to inefficient use of the network resources 
0. Therefore, there has been increasing interest in designing better client association policies 
BU— Q. Whereas the previous approaches are hard to apply in 60 GHz wireless access networks 
due to the special characteristics of the 60 GHz channel, and the differences with the rest wireless 
access technologies ff8l--|fTT| (namely, severe channel attenuations, high path loss, directionality, 
and blockage), novel mechanisms must be designed to provide optimal resource allocation. Our 
previous approach [fl2~l was the first to study the client association in 60 GHz wireless access 
networks. However, the focus was on the network performance (achieving load balancing) and 
not on optimizing the benefit of the individual clients. 

This paper considers the special characteristics of the 60 GHz access channel and poses the 
client association optimization problem, where the objective is to maximize the total clients 
benefit in the network. To address the problem, we propose an iterative approach that combines 
two auction algorithms. We compare our solution method to basic association policies, already 
in use in the present 60 GHz communication technologies under standardization (802.15.3c, 
802. Had) 031. 

The rest of the paper is organized as follows. A description of the system model and the 
problem formulation is presented in § HB In § HEH we describe the solution approach to the 
multi-assignment problem. In § [IV] numerical results are presented. Lastly, § |V] concludes the 
paper. 

II. System Model and Problem Formulation 
We consider a mmW network where m access points (APs) that can serve n clients and 
n > m. An AP i can serve more than one client. Moreover, every client j must be associated to 
just one AP. The set of clients to which AP i can be assigned is a nonempty set A(i). Moreover, 
we introduce the set B(j) as the nonempty set of APs that can serve client j. An assignment S 
is defined as a set of AP-client pairs (i,j), with j E A(i), where each AP i can be part of more 
than one pair (i,j) E S, and where every client j must be part of only one pair (i,j) E S. An 




Figure 1: Example mmW wireless access network. 

illustrative example of access network is shown in Figure [Q where the clients positioned inside 
a disc with radius r centered at the location of AP i can be associated with that AP. 

Every node is equipped with steerable directional antennas and it can direct its beams to 
transmit or to receive fl8]. We assume that AP i can support its clients with a separate transmit 
beam. We consider the case where all receiver nodes are using single-user detection (i.e., a 
receiver decodes each of its intended signals by treating all other interfering signals as noise) 
and assume that the achievable rate from AP i to client j G A(i) is 

where W is the system bandwidth, P^ is the transmission power of AP % to client j, Gij is 
the power gain from AP i to client j, N is the power spectral density of the noise at each 
receiver, and Ij is the interference spectral density at client j. We use the Friis transmission 
equation together with the flat-top transmit/recieve antenna gain model [9|, where a fixed gain 
is considered within the beamwidth and zero gain is considered outside the beamwidth of the 
antenna. We capitalize on the well studied 60 GHz propagation characteristics [9|, such as 
highly directional transmissions with very narrow beamwidths and increased path losses due to 
the oxygen absorption, in order to assume that the communication interference Ij is very small 
and does not affect significantly the achievable communication rates in the network. We remark 
that all the assumptions that we make above are natural for 60 GHz ||9l. 

We denote by Qj the demanded data rate of client j. The benefit of client client j that 
is associated with AP i is given by the ratio of Rij/Qj. The general objective is to find an 
assignment that maximizes the sum of such benefits, namely the total benefit of the network. 
Therefore, the association problem is modelled by the following linear optimization problem 



max 22 ~n~ Xi i ( 2a ) 
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s.t. 



^2 Xi i - 1 ' Vi = 1, . . . , m, (2b) 



jeA(i) 



J^ x tj = 1, Vj = l,...,n, (2c) 

%>0, V(i,j)eC (2d) 

The objective function of © is the total network benefit, where C is the set of all possible AP- 
client assignment pairs (i,j) (note that S is a subset of C) and (:%)jeA(i) are binary decision 
variables, indicating the client association. In particular, Xij = 1 if client j is associated to AP 
i and Xij = otherwise, for all i and j £ A(i). (l2bl) and (|2cl) ensure that each AP will be 
assigned to one or more clients and each client will be associated to one AP. Note that from 
an assignment S, we can potentially recover a solution to problem (J2j by setting Xij = 1 if 
(i,j) € S and x i: j = otherwise. An assignment that gives a feasible solution to problem © 
is therefore defined as feasible assignemnt. In what follows, we present the proposed solution 
approach. 

III. Solution Approach 

The considered problem (0) is a classical multi-assignment problem, where an AP can be 
assigned to more than one client. Unfortunately, there are no specialized network flow methods 
that can efficiently solve this class of assignment-like problems. There are approaches that apply 
general purpose network methods such as primal-simplex, primal-dual, or relaxation methods, 
which may have high complexity [[14|. Moreover, general methods for linear optimization, 
such as the simplex or even interior point methods, do not exploit the particular structure 
of the considered multi-assigment problem at hand (a network optimization structure) and are 
not amenable for distributed computation. Thus they are generally less efficient than network 
optimization methods [[131 . Consequently, we resort to network optimization theory and propose 
a solution method that combines auction algorithms to solve efficiently problem (0). 

We start by converting problem (0) into a typical minimum cost flow problem [[151 by 



introducing a virtual supernode s that is connected to each AP 
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min ^ ~FT~ Xi J ( 3a ) 
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s.t. 2_. x ij ~ x si = 1) Vi = 1, . . . , m, (3b) 



jGA(i) 



^ x si = n - m, (3c) 



i=l 



^2 x ij = !> Vj = l,...,n, (3d) 

ies(j) 

z<i>0, V(i,j)eC, (3e) 

a^si > 0, Vz = 1, . . . ,n (3f) 

where the sign of the benefit was reversed (cost coefficient) compared to problem d2), minimiza- 
tion replaced the maximization and x^ was extended to include also the supernode s. By using 
the terminology of network optimization, x^j has the meaning of amount of flow between i and 
j, and the first constraint ensures that the flow supply of each AP i is one unit, while the second 
one declares that s is the source node and the flow that generates is of n — m units. Therefore, a 
flow of one unit will reach each client j. The last two constraints declare that the flow of each 
arc may be infinite, where an arc between i and j denotes the connection (i,j). A solution to 
the minimum cost flow problem © is the same to the initial multi-assignment problem ©. 

By using the duality theory for minimum cost network flow problems [fTSl §4.2] we formulate 
the dual problem 



m n 



min \ 7Tj + \ pj + (n — m)X (4a) 

aXn+Pj >fy V(i,j)eC, (4b) 

A > 7Tj, Vz = l,...,m (4c) 



where — 7Tj is the Lagrangian multiplier associated with constraint (|3bl representing the price of 
each AP i, X is the Lagrangian multiplier associated with constraint (l3cl) representing the price 



'We consider a network where supernode s generates n — m units of traffic and is connected to each AP j by a zero cost 
arc (s,i). The traffic that is generated at each AP (supply) is of one unit. AP i is connected to client j by an arc (i,j) with 
cost —Rij/Qj. 



of the supernode s (recall that s is the source of the flows), and pj is the Lagrangian multiplier 
associated with constraint (l3dl representing the price of each client j. The optimal solution to 
problem © allows us to derive the optimal solution to © [15, §4.2, §5]. 

In order to solve problem © we need some technical intermediate results. We start by giving 
the definition of e-Complementary Slackness (e — CS): Let e be a positive scalar, we say that 
an assignment S and a pair (ir,p) satisfy e — CS if 

Ki+Pj >%*--*, V(i,j)eC, (5a) 

*i+Pi = ^-, V( 8) j)65, (5b) 

7Tj = max 7r fc , Vi s.t. i has more than one pair (i,j) G 5 1 (5c) 

fc=l,...,m 

Proposition 1: Consider problems © and ©. Let 5 be a feasible solution for problem © 
and consider a dual variable pair (n,p). Let e < 1/m and assume Rij/Qj be integer Vz,j. If 
e — CS 1 conditions © are satisfied by S and n,p, then S 1 is optimal for problem ©. 

Proof: The proof is ad-absurdum. Assuming that S is not optimal, then there is a new 
assignment that can improve the objective function © and can give us a new solution: Let E 
be a cycle, namely a collection of arcs that start and end with the same node, that includes also the 
supernode s: E = (s, h,J2, ii-, ■■■■, ik-i,jk, ik, s). In this solution, the nodes i t represent the APs, 
while the nodes j t represent the clients and (it,jt) G S,j t G A(i t _i), (it-i,jt) £ S,t = 2, ..., k. 
Based on max-flow theory [fT5l §5], augmentation along E is achieved by replacing (it,jt) G S 
by (it-i,jt) in S, t = 2,...,k. AP ik must be assigned to more than one clients prior to the 
previous operation because the arc (ik,jk) will exit the assignment and therefore, the AP i k will 
be left unassigned. This will result to an infeasible solution to problem ©. Moreover, k < m 
since E cannot contain repeated clients. Considering also that e < 1/m we conclude that ke < 1. 

Since we achieved strict cost improvement in the previous operation, we have 

E^ + ^Ef^ < 6 > 

t=2 ^'* t=2 ^ 3t 

In order to reveal the e — CS 1 conditions ©, we transform © as 

k ' R \ - fc , / f? 



E 3*-** + ><E ^-ft. • (?) 



t=2 v ^ Jt 7 i=2 v ^-?* 

Now using the e — CS conditions ©, © can be written as 



Algorithm 1 Forward Auction for Client Assignment Algorithm 2 Reverse Auction for Client Assignment 

Require: S, (%,p) and A from forward auction 
Ensure: (1) n + pj > Rij/Qj - e, V(i, j) e C and 
(2) m +pj = Rij/Qj, V(i, j) 6 S 
while there are unassociated clients do 
client j is unassociated in S 



Require: Initial values of S,p 

Ensure: R^/Qj -pj > max teA(l) {Rij/Qj - p k } - 
e, V(u]£S 

while there are unassigned clients do 
client j is unassigned in S 
find the best client ji such that: 
j % = argmax, e A(j) {Rij/Q - Pj} , 
m = max j£ A(i) {Rij/Qj ~Pj}, 



v j'e 



A{i),j^ji { R v/Q] ~Pj}: 



if ji is the only client in A(i) then 

LUj — > — OO 

end if 

hji = Ph + 1M -uJi+e = Rij/Qji -cJi + e 
Pj = maxj e p(j-) bij, where P(j) is the set of APs 
that client j received a bid, 
remove any pair (i,j), where j was initially 
assigned to some i under S, and add the pair (ij,j) 
to S with ij = argmaxjgpQ-) fri- 
end while 



find the best AP ij such that: 
ij = argmaxjgB^) {Rij/Qj - n} , 
Pj = maxjgsQ) {Rij/Qj - 7r 4 } , 
ujj = maxizBW&ij {Rij/Qj ~ Ki} , 
if ij is the only AP in B(j) then 

LOj — > — OO 

end if 

S = min {A — 7Tj- , /3j — Wj + e} 

add (ij, j) to S: pj = ^ - <5, tt^ = tt.^. + 5 

if 6 > then 

remove the pair (ij t j id) where j id was ini- 
tially assigned to ij under 5 

end if 
end while 



k 



7T; 



-tti t j t 



Qj 



fe-i 



i<E ^-^ME^+(*-i) e . 



(8) 



t=2 i=2 v ^• ?t ' t=l 

From ([8]) we have 1 — (k — l)e < 7Tj 1 — 7r ifc which contradicts fee < 1, because AP ik is assigned to 
more than one clients, i.e., n ik > n^, [compare with (|5cT )1. We conclude that our first assumption 
on that S is non optimal is wrong, which implies that S is optimal. We can get similar results 
considering that supernode s is not part of E. ■ 

Note that in general, Rij/Qj is not an integer as required by the Proposition [T] i.e., rounding 
those to the closest integer value or scalled to an integer value is necessary before running 
the algorithms. However, in 60 GHz access networks, the effect of rounding influences slightly 
the true optimal value of Problem dU, because we can assume Rij ^> Qj, and as a result the 
fractional part of Rij/Qj is relatively smaller than the integer part of Rij/Qj. 

Based on Proposition [fl we are now in the position to present the solution method to prob- 
lem ((4]) by an auction mechanism. First, a forward auction algorithm associates each AP to one 
client, see Algorithms [Q Then, a modified reverse auction is applied to assign the rest of the 
clients to the available APs, see Algorithms [Q Finally, we show that the execution of the two 
algorithms terminates with an optimal solution by a finite number of iterations. 

In particular, we start from a feasible assignment S and the corresponding (ir,p) pair that 
satisfy the first two e — CS conditions. We apply Algorithm \T\ until each AP is associated 
with a single client and until the e — CS conditions are satisfied. At this stage, some of the 



clients can still be unassigned. We then apply Algorithm |2] that gets as input the assignment 
achieved by Algorithm [TJ (S and (%,p)). We compute the maximum initial profit for the APs 
A = max i=1 m 7Ti. The iterative Algorithm |2] maintains an assignment S, where each AP is 
associated with at least one client, and a pair (ir,p) that satisfies the first two e — CS conditions. 
Algorithm |2] terminates when all unassigned clients have been assigned to an AP. While A is 
kept constant through the execution of Algorithm |2] (l5cT > is satisfied upon termination. 

Proposition 2: Consider Algorithm [TJ and El Let Algorithm [TJrun first and then let Algorithmic] 
run iteratively. Algorithm |2] terminates in a finite number of iterations with an optimal AP-client 
assignment when e < 1/m. 

Proof: In order to prove the optimality and the convergence of the modified reverse auction 
algorithm we have to show that a) The modified reverse action algorithm iterates by satisfying 
e — CS conditions © and A = max i= i ... m tTj, b) The algorithms terminates after a finite number 
of iterations with a feasible assignment. 

The proof of a) is a straight forward application of the well known theory for auction 
algorithms [fT5l §7] to show that if the e — CS conditions © and A = max i= i v .., m 7r; are 
satisfied at a start of an iteration, they are also satisfied at the end of the iteration. 

To show b), we observe that an AP i can receive a bid only a finite number of times after 
7Tj = A. This is true due to that in each iteration the corresponding client will be assigned to 
AP i without changing the association of already assigned clients to AP i (see Algorithm |2]>. 
Lastly, at the end of each iteration when AP i receives a bid, the profit 7Tj is either equal to A or 
else increases by at least e. Since A is an upper bound in the profits throughout the algorithm, 
the main outcome is that each AP can receive a finite number of bids (finite termination). ■ 

Proposition [2] is very important to get an insight of the behavior of the proposed algorithm 
in general networks, see Section [IV] for numerical examples. 

IV. Numerical Analysis 
In this section we present a numerical evaluation of the proposed algorithm in a multi-user 

multi-cell environment. We compare the proposed solution approach to a) random associa- 
tion policy and b) RSSI-based policy, which is the standard association mechanism used in 
802. 11 ad [[H networks. 

We define the SNR operating point at a distance d form any AP as 




Number of clients 




Figure 2: Total benefit © vs. number of 
clients (10 APs). 



Figure 3: Total benefit © vs. number of 
APS (100 clients). 



SNR(d) 



P Q \ 2 /(lQiT 2 N W) d<d 

P \ 2 /(16tt 2 N W) ■ (d/d y v otherwise . 
We consider circular cells (assumed for simplicity, without loss of generality), as depicted 
in Figured] The radius of each cell r is chosen such that SNR(r) = 10 dB. APs are located such 
that the distance between any consecutive AP is D — l.lr. The clients are uniformly distributed 
among the cells, and the potential AP-client association is found as pointed out in Figure [Q 

We set A = 5 mm, iVo = -134dBm/MHz, W = 1200 MHz, and d = lm, see ©. Moreover, 
we set Ij = 0, Pij = Pq = 0.1 mW and Gjp = G^ x = 1. Furthermore, we assume that QjS are 
uniformly distributed on [0, 100] Mbits/s. The algorithms were implemented in MATLAB and 
run on an Intel Core 2 Duo 2.40 GHz processor with 8 GB RAM. 

Figure |2] depicts the total benefit ^ R^jQj (main objective in problem ©) achieved 
by our solution approach in comparison to the optimal solution of ©, the received signal 
strength (RSSI) based mechanism (adopted by 60 GHz standards) and the random association 
methodology, where 10 APs are present in the network and the number of the supported clients 
varies. We observe that the auction-based approach achieves optimal performance and improves 
the performance of RSSI-based mechanism up to 75% (especially in high load conditions). 
Figure |3] depicts the total benefit in the network when the number of APs varies for fixed 100 
clients. The behavior of our approach is similar to Figure |2] evincing its optimal and scalable 
performance. 

Figure @] shows the performance of the proposed approach for different network sizes (20 
clients 2 APs, 40 clients 4 APs, etc) and parameters e. Note that we have e < 1/m for all 
considered cases, in order to guarantee optimal performance, see Proposition [Q Results show 
that the convergence time of the proposed association algorithm is approximately linear in n. 
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Figure 4: Time (sec) vs. size of the network. 

V. Conclusion 
We considered the problem of optimizing the allocation of the clients to APs in mmW wireless 

access networks. The objective in our problem formulation was to maximize the total clients 

benefit. We presented a solution approach based on forward and reverse auction algorithms. 

Both theoretical and numerical results evinced the optimal, scalable and time efficient behaviour 

of our approach. Thus, it could be well applied in the forthcoming 60 GHz wireless access 

networks. 
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